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ABSTRACT 

We present a new estimate of the lithium abundance in the atmosphere of the brown 
dwarf LP 944—20. Our analysis is based on a self-consistent analysis of low, intermedi- 
' ate and high resolution optical and near-infrared spectra. We obtain log N(Li) = 3.25 

\^ , ± 0.25 using fits of our synthetic spectra to the Li I resonance line doublet profiles 

■ observed with VLT/UVES and AAT/SPIRAL. This lithium abundance is over two 

orders of magnitude larger than previous estimates in the literature. In order to obtain 
good fits of the resonance lines of K I and Rb I and better fits to the TiO molecular 
absorption around the Li I resonance line, we invoke a semi-empirical model atmo- 
\ sphere with the dusty clouds located above the photosphere. The lithium abundance, 

however, is not changed by the effects of the dusty clouds. We discuss the implications 
of our estimate of the lithium abundance in LP 944— 20 for the understanding of the 
\ properties of this benchmark brown dwarf. 

^ . Key words: stars: individual: LP 944— 20 - lithium abundance - stars: fundamental 

parameters - stars: late-type - stars: brown dwarfs - stars: evolution 
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1 INTRODUCTION 

M dwarfs are of special interest to many branches of modern 
astrophysics. The population of these numerous low-mass 
stars (M ^ 0.6 M0 ) can contain an appreciable amount of 
the baryonic matter in the Galaxy. The coolest M dwarfs 
(later than M6) could be young brown dwarfs (Martm , Re- 
bolo & Zapatero Osorio 1996) if they are young enough. The 
existence of brown dwarfs was predicted by Kumar (1963) 
and Hayashi & Nakano (1963). The first unambiguous brown 
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072.C-0110(B), the Anglo-Australian Telescope at Siding Springs 
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strument; and the Keck Observatory in Mauna Kea, Hawaii. 



dwarfs Teide 1 and Gliese 229B were reported by Rebolo, Za- 
patero Osorio & Martm (1995) and Nakajima et al. (1995), 
respectively. 

The "Lithium test" observation of the Li I resonance 
doublet at A 670.8 nm was proposed by Rebolo et al. (1992) 
and developed by Magazzii et al. (1993) to distinguish young 
brown dwarfs from very low-mass stars. This test has allowed 
the identification of many brown dwarfs of spectral classes 
M and L (Rebolo et al. 1996; Martm et al. 1997, 1999; Kirk- 
patrick et al. 1999). 

Due to the low temperatures in the cores of brown 
dwarfs of masses M < 55Mjup lithium is not destroyed. 
From 55Mjup < M < 75Mjup, hthium is destroyed on 
timescales that can be used to obtain nuclear ages for these 
very low-mass objects (e.g. Chabrier et al. 2000). However, 
the application of the lithium test requires improvement of 
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our knowledge about the physics of atmospheres in low tem- 
perature regimes (Pavlenko et al. 1995). 

LP 944-20 (other names are APMPM J0340-3526, BRI 
B0337-3535, LEHPM 3451, 2MASS WJ0339352-352544) is 
an archetypical brown dwarf of spectral type M9. A dim 
red dwarf-hke star, Luyten Palomar (LP) 944-20 (M9 V) 
was cataloged by Luyten & Kowal (1975), and the SIMBAD 
database lists more than 125 references for this object up to 
May 2007. It is interesting to note that it is comparatively 
bright (Mboi = 14.22, Dahn at al. 2002) and also appears to 
be a young object, a fact that is confirmed by the lithium 
detection of Tinney (1998). 

Using evolutionary tracks Tinney (1998) estimated 
the age of LP 944-20 to be between 475 and 650 Myr. 
Ribas (2003) found some evidence for the membership of 
LP 944—20 in the Castor moving group in the solar vicinity 
(J? ~ 5 - 20 pc from the Sun). The age of LP 944-20 as a 
member of the Castor group LP922-20 has been estimated 
as 320 ± 80 Myr (Ribas 2003). 

An unexpectedly strong, solar like X-ray flare (log 
{Lx/Lhoi = -4.1) of duration ~ 2 hours was registered by 
Chandra (Rutledge et al. 2000), but later XMM-Newton ob- 
servations by Martin & Bouy (2002) did not detect any X- 
ray activity. The quiescent X-ray luminosity of LP 944— 20 
is lower than that of the active Sun, and X-ray flares 
on LP 944—20 are very rare events in comparison with 
more massive M-dwarfs (Hambaryan et al. 2004). How- 
ever, LP 944—20 is not the only M9 dwarf showing X-ray 
flares. Recently Stelzer (2005) observed an X-ray flare (log 
{Lx/Lhoi = -3.3) from the brown dwarf binary with spectral 
types M8.5/M9 G1569 B, which has an age of about 300 Myr 
(Zapatero Osorio et al. 2004, see also Viti & Jones 1999). 

Our knowledge about the structure of the outermost 
layers of ultracool dwarfs is far from complete (e.g. see dis- 
cussion in Liebert et al. 2003). Indeed, Berger et al. (2001) 
reported the discovery of both quiescent and flaring radio 
emission from LP 944—20. They observed radio luminosities 
of 80 /iJy in quiescence and 2 /xJy at flare peak respectively, 
which are several orders of magnitude larger than those pre- 
dicted by the empirical relations between X-ray and radio 
luminosities that has been found for late-type stars. From 
the comparison of the continuous and flare fluxes at 8.5 GHz, 
Berger et al.(2001) concluded that the dwarf has a rather 
weak magnetic field in comparison with fiaring M-stars, and 
its emission mechanism is the most likely non-thermal syn- 
chrotron radiation. 

Tsuji et al. (1996) predicted "dusty effects" in the atmo- 
spheres of cool dwarfs with Tcff < 2600 K. The presence of 
large amounts of dust in the highly dynamical atmospheres 
of fast rotating ultracool dwarfs can yield "weather phenom- 
ena" (Jones & Tsuji 1997, Gelino et al. 2001; Martin, Zap- 
atero Osorio & Lehto 2001; Clarke et al. 2002; Bailer- Jones 
& Lamm 2003, Caballero et al. 2003; Koen et al. 2004). 
LP 944—20 also shows significant rotational broadening in 
high-resolution spectra. The v sini values determined by dif- 
ferent authors are in the range of 30 ± 2 km/s (Tinney & 
Reid 1998; Mohanty & Basri 2003, Guenther & Wuchterl 
2003; Jones et al 2005; Zapatero Osorio et al. 2006). The ro- 
tation period has not been determined yet, but photometric 
variability was reported by Tinney & ToUey (1999). These 
surface inhomogeneities could be due to magnetic spots or 
dusty clouds modulated by stellar rotation. 



Tinney (1998) reported the observation of a strong 
lithium line in the spectrum of LP 944—20 as evidence of its 
youth and substellar nature, and estimated a lithium abun- 
dance of log N(Li) = 0.0 ± 0.5, indicating it had depleted 
an important amount of lithium. This is so far the only de- 
termination of the lithium abundance in LP 944—20, which 
was obtained via comparison of the observed feature width 
with theoretical curves of growth. However, the observed 
atomic line at 670.8 nm is measured with respect to the 
molecular background formed by TiO lines and therefore a 
detailed spectral analysis needs to be made. Here we present 
new computations using different approaches and model at- 
mospheres indicating that LP 944—20 has not depleted any 
lithium. 

This paper is organized as follows: Section 2 presents the 
spectroscopic data of LP 944—20 available in the literature. 
Section 3 describes the procedure followed to compute syn- 
thetic spectra and spectral energy distributions using stan- 
dard model atmospheres. We fit these theoretical spectra to 
the observations. In Section 4 we discuss the results obtained 
in the framework of the conventional approach and we de- 
rive a new estimate of the lithium abundance in LP 944—20. 
In section 5 we investigate the impact of the possible de- 
viations of physical parameters on our results. We propose 
a semi-empirical (SE) model for the formation of molecular 
features and strong resonance lines of alkali elements in the 
spectrum of LP 944—20. We present the results of applying 
the SE model for fitting K I, Rb I and Li I lines in high- 
resolution spectra and spectral energy distributions. We de- 
rive a new estimate of the lithium abundance in LP 944—20 
with our SE model atmosphere and in Section 6 we discuss 
the implications of our results. 

Coloured versions of plots are available at 
ftp://star-ftp.herts.ac. u k/pub/Pavlenko/lp944.20] 



2 OBSERVATIONS 

The observed spectra used for this paper have been pub- 
lished elsewhere. In order to simplify the text of the paper 
we adopt a labelling syntax for them as follows: 

• CASPEC: optical spectra obtained by Tinney & Reid 
(1998) using the Cassegrain Echelle Spectrograph on the 
ESO 3.6 m telescope (R ~ 18000). The infrared part of the 
spectral energy distribution was obtained by UKIRT (see 
references in Leggett et al. 2001). The combined spectrum 
was taken from the Leggett 's database 

(ftp:/ /ftp.jach.hawaii.edu/pub/ukirt/skl/ ). 

We believe Leggett et al. (2001) normalised the spectral data 

using broad-band photometry but the error on this is of the 

order of 10%. We have rescaled the flux beyond 1 ^m by 

10% to improve the match between observed and synthetic 

spectra. 

• SPIRAL: spectra of medium resolution obtained with 
SPIRAL Phase A fibre-fed spectrograph on the AAT. The 
FWHM was typically 2.4 to 2.7 pixels over the whole data, 
the GFWHM, i.e. FWHM of the smoothing Gaussian, was 
of order 0.883 A for the Li I 670.8 nm line (R=7500). See 
Kenworthy et al. (2001) for the data reduction details. 

• NIRSPEC: the near-infrared spectrum of LP 944-20 
(1.240-1.258 ^im , R ~ 22000) was taken with Keck II. AU 
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details about those observations have been pubhshed in Za- 
patero Osorio et al. (2006), and Martin et al. (2006). 

• UVES: high spectral resolution spectrum obtained with 
UVES (Ultraviolet- Visual Echelle Spectrograph) on the 
VLT Unit telescope 2 (KUEYEN) at Paranal in service 
mode. The setting simultaneously covers the wavelength re- 
gions from 667.0 to 854.5 nm. We used the average of 15 
individual spectra obtained for the purpose of radial veloc- 
ity monitoring. Details of the data reduction and calibration 
can be found in Guenther & Wuchterl (2003). The theoreti- 
cal resolution of our UVES spectra is of order 40000, corre- 
sponding to the GFWHM = 0.17 A. 



3 STANDARD SPECTRAL MODELLING 

To determine the surface temperature (Tcff) of LP 944— 20 
we have used the equation L = iTvR^aT^g, normalized to so- 
lar units. The trigonometric parallax of LP 944—20 is found 
to be 200.7" ± 4.2" , which is the mean value of the measure- 
ments obtained by Tinney (1998) and Dahn et al. (2002), 
yielding a distance of 5.0 pc to the dwarf. By using the 
bolometric corrections for the J- and /("-bands and the ob- 
served near-infrared photometry of LP 944—20 available in 
the literature, we have derived its luminosity at log L/Lq = 
—3.79 ± 0.03, which is in agreement with the previous deter- 
mination by Dahn et al. (2002). The error bar accounts for 
the uncertainties in the photometry, bolometric corrections 
and parallax. 

The radii of "old" (S5500 Myr) brown dwarfs are nearly 
independent of mass and age, with a mean radius of 0.09 Rq 
(see Burgasser 2001) we derive Tcft = 2170 K. However, 
the fact that lithium is indeed detected in the optical spec- 
trum of LP 944—20 indicates that this object is significantly 
younger than other field dwarfs of similar classification, 
which do not show the lithium feature in absorption. Tinney 
(1998) and Ribas (2003) state that the age of LP 944-20 is 
in the range 240-650 Myr, with a likely value at 320 Myr. For 
such a young age, recent evolutionary models predict brown 
dwarf (0.05-0.07 Mg) radu of about 0.10 R@ (Chabrier & 
Baraffe 2000; Burrows et al. 1997), i.e. 11% larger than typ- 
ical field brown dwarfs. This yields a temperature Toft = 
2040 K, which is the value that we will use throughout the 
present paper. Recently, Bihain et al. (2006) have shown that 
proper motion Pleiades brown dwarfs may have radii simi- 
lar to those of much older substellar objects, suggesting that 
brown dwarfs collapse faster than expected. To account for 
this, we estimate the uncertainty in our TcS determination 
of LP 944— 20 to be ±150 K, which includes the luminos- 
ity error bar and the various object sizes valid for an age 
interval 120-1000 Myr. 

Our strategy is divided into the following steps; (i) The- 
oretical spectra are computed and fitted to the observed 
spectral energy distribution of LP 944—20 from 0.65 to 2.5 
nm . (a) We then compare our synthetic spectra obtained 
following a standard classical approach with some atomic 
and molecular features in the observed spectrum (in partic- 
ular, the resonance lines of K I, Rb I) . Here we use spectra 
of medium and high spectral resolution. Finally, (in) we de- 
termine the sensitivity of our best fit lithium abundance to 
the input parameters of our models. 
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Figure 1. Molecular densities of lithium containing species in 
the 2000/4.5/0.0 model atmosphere {T^a = 2000 K, log g = 4.5, 
solar metallicity). 



3.1 Synthetic spectra 

Theoretical spectral energy distributions (hereafter we use 
the term "synthetic spectra") were computed for model 
atmospheres of dwarfs with effective temperatures T^g = 
1800-2400 K and log g = 4.0 - 4.5 (cm s"^) from the 
DUSTY and COND model atmospheres (AUard et al. 2001) 
of solar metallicity (Anders & Grevesse 1989). Hereafter we 
use the model atmosphere notation 2400/4.5/0.0 to mean 
Tcff = 2400 K, log g = 4.5, [M/H] = 0. We wiU focus on 
results obtained with 2000/4.5/0 model atmosphere because 
as previously indicated, 2040 K is an appropriate effective 
temperature of LP 944—20, and log g = 4.5 since LP 944-20 
is considered to be younger than field dwarfs of similar clas- 
sification. Nevertheless, computed spectra for other temper- 
atures and gravities were used to study the dependence of 
our results on these parameters. 

The dominant opacity sources in the optical and in- 
frared spectra of LP 944—20 are absorption by band systems 
of diatomic molecules, such as TiO and VO. Computations 
of synthetic spectra were carried out by the program WITA5 
(Pavlenko et al. 2000) assuming LTE, hydrostatic equilib- 
rium for a one-dimensional model atmosphere and without 
sources and sinks of energy. The equations of ionisation- 
dissociation equilibrium were solved for media consisting of 
atoms, ions and molecules. We took into account ~ 100 com- 
ponents (Pavlenko 1998a). The constants for the equations 
of chemical balance were taken mainly from Tsuji (1973) 
and Gurvitz et al. (1989). 

Our lithium containing species list includes Li I, Li II, 
LiOH, LiH, LiF, LiBr, LiCl. Our computations show that 
within our {Tos , log g) range neutral lithium dominates 
across the atmospheres - only 5% of lithium atoms are bound 
in molecules in the outermost layers of atmosphere (Fig. [l|. 
The most abundant lithium containing species are molecules 
of LiCl (in the outermost layers) as well as Li II and LiOH 
(in the photospheric layers). 
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Figure 2. A plot identifying the main features in a 2400/4.5/0.0 
model spectrum, showing their relative strengths. 



3.2 Molecular and atomic lines opacity 

The molecular line data are taken from different sources: 
TiO line lists by Plez (1998); CN lines from CDROM 18 (Ku- 
rucz 1993); CrH and FeH lines from Burrows et al. (2002a) 
and Dulick et al.(2003), respectively. Atomic line list data 
are taken from VALD (Kupka et al. 1999). 

The profiles of molecular and atomic lines are deter- 
mined using the Voigt function H{a, v) except for the strong 
resonance doublet lines of Na I (0.5891, 0.5897 /^m ) and K I 
(0.7667, 0.7701 ^im ). (see the section [3]4} . The parameters 
of natural broadening C2 and van der Waals broadening Cg 
of absorption lines are taken from Kupka et al. (1999), or 
in their absence computed following Unsold (1955). In our 
computations we used a correction factor E for van der Wals 
broadening of atomic lines. Parameter E describes differ- 
ences in broadening parameters computed in the framework 
of classical and quantum physics approaches. For resonance 
lines of alkali metals E — 1-2 (see Andretta et al. 1991). Ow- 
ing to the low electron and ion densities in low temperature 
M dwarf atmospheres. Stark broadening may be neglected. 
In general the effects of pressure broadening prevail. Compu- 
tations for synthetic spectra to be fitted to observed spectra 
across the 0.65 - 0.9 fim are carried out at intervals of 0.5 A. 
For comparison with high resolution spectra we computed 
spectra with an interval of 0.02 A. The spectrum broadening 
is modelled by Gaussian plus rotational profiles set to the 
resolution of the observed spectra. Rotational broadening 
was taken into account following the Gray (1976) formulae 
with V sini = 30 km/s (see Jones et al. 2005). VO band opac- 
ity is computed in the frame of the JOLA approximation (see 
Pavlenko et al. 2000 for more details). The relative impor- 
tance of the different opacities contributing to our synthetic 
spectra is shown in Fig. [5] 

In this paper we also pay careful attention to the mod- 
elling of the rubidium resonance line at 780.2 nm. Resonance 
lines of Rb I have superfine structure (Lambert & Mallia 
1968, Lambert & Luck 1976), which is important to account 
when considering weak absorption lines, but the Rb I lines 
in our spectra are strong. Moreover, they are broadened by 
rotation and turbulent motions and so the superfine split- 
ting of Rb I lines can be neglected (see ibid.). We assume 



that all atoms of rubidium are the ^*Rb isotope. Most fea- 
tures of the observed spectra should be fitted using the same 
atmospheric parameters of temperature and gravity. We will 
use our fits to the Rb I line as a test of the consistency of 
our spectroscopically derived results to log g and lithium 
abundance in the atmosphere of LP 944—20. 



3.3 Observed profiles of Li I and Rb I 

The fiux variations across Li I and Rb I line regions obtained 
in a few separate observing runs are shown in Fig. |31 We 
show there a sequence of the first 5 observed fluxes from a set 
of 15 spectra which were used to get the combined spectrum 
used in our analysis. It is worth noting that changes in the 
Rb I line 780.2 nm profile are rather small. However, we 
see here some variations in molecular background formed 
by the haze of TiO lines. In general, the response of the 
TiO molecular densities to changing temperatures in the 
line forming region should be much higher in comparison 
with changes of the ionisation equilibrium of alkali metals 
(see Pavlenko 1998). Some variations of the observed fluxes 
across the 780 nm region are real. Analysis of the topic is 
beyond the scope of our paper. 

For the 670 nm region, we find that the data from indi- 
vidual UVES epochs appear to be reasonably fit, albeit at 
low signal-to-noise, by the same parameters as the higher 
signal-to-noise Rb lines. However, when the dataset is com- 
bined (15 low signal-to-noise epochs spread over more than 
a year) a considerably broader Li I profile results. This is 
not due to radial velocity changes for which this dataset has 
already being corrected to a level of the order of 80 ms~^ 
(Guenther & Wuchterl 2003). This apparently changing line 
profile is presumably due to either variability and/or rela- 
tively low signal-to-noise ratio of the UVES data at these 
wavelengths and leads us to fit the combined profile with a 
lower effective resolution. We find the lithium region to be 
best fit with v sini = 32 km/s and a GFWHM = 1 A. 



3.4 K and Na resonance lines 

In the spectra of ultracool dwarfs the resonance lines of Na 
I and K I are very strong. They govern the spectral energy 
distributions of L-dwarfs across a wide spectral region (see 
Pavlenko 2001 for more details). Their formally computed 
equivalent widths may be of the order of a few thousand 
A. For these features we cannot use a classical collisional 
approach to compute the wings of these superstrong lines. 
In the dense, cool atmospheres of late M, L and T dwarfs 
the pressure broadening of K I and Na I can be computed 
with a quantum chemical approach (Burrows & Volobuyev 
2003; AUard at al. 2003). 

For this work we use potentials of quasi-stationary 
chemical interactions of atoms K and Na with the most nu- 
merous species, i.e., atoms He and molecules H2 computed 
by GAMESS (Granovsky et al. 1999) . Our procedure is de- 
scribed in more detail in Pavlenko et al. (2007). Here we 
point out that in computations of K I profiles we used a 
combined profile: cores of these lines were computed in the 
frame of the collisional approach, and their wings {SX > 40 
A) were treated by quasi-stationary theory. It is worth not- 
ing that K and Na lines have comparatively weak wings in 
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our range of Tea and log g. Thus, in observed spectra we see 
only cores of strong absorption lines against a background 
of molecular bands. 

For other alkali lines we used a conventional, i.e., col- 
lisional theory of pressure broadening, using van der Waals 
formulae for computation of damping constants described 
in Pavlenko (2001). The relative strength of the resonance 
lines of neutral alkali lines is shown in Fig. [2] 

3.5 Fit to the observed spectral energy 
distribution 

To get the best fits of our theoretical spectra to the ob- 
served spectral energy distribution we followed the scheme 
described in Pavlenko et al. (2006) and references therein. 
Namely, we find the minimum value of the following equa- 
tion: 

S{h, h) = E(/h X H^y-' - H'"°n' (1) 

Here H°^^ and 7^^°* are observed and computed fluxes, /h 
is a normalisation parameter, fs is the relative wavelength 
shift of the observed spectra, and fg is related to the broad- 
ening parameters. In our case, we applied both rotational 
and instrumental broadening. These were modelled by Gray 
(1996) formulas (rotational broadening) and Gaussian func- 
tions (instrumental broadening). 



4 RESULTS BASED ON STANDARD 
SPECTRAL MODELLING 

4.1 Fit to the observed spectral energy 
distributions 

In Fig. ID we show fits to the observed spectral energy distri- 
butions (SEDs) of LP994-20 across the 0.65-2.5 ^m region. 
We carried out our fits for the DUSTY and COND model 
atmospheres. We exclude from the fitting the spectral region 
0.86-1.96 fj.m because there is a gap in the observed spec- 
trum, and strong spectral features formed by absorption of 
H2O and other molecules at these wavelengths cannot be 
fitted properly even for the objects of higher Tcft (Pavlenko 
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Figure 4. Fits of the synthetic spectra computed for DUSTY and 
COND model atmospheres 2000/4.5/0.0 to the observed spectral 
energy distribution of LP 944—20. There is no data between 0.9 
and 1.1 ^m. 



et al. 2006). Nonetheless, our approach allows us to use the 
optical region governed by TiO bands where flux is very sen- 
sitive to Toff , and the region beyond 1.96 /im dominated by 
H2O and CO bands. Line lists of these molecules for the 
latter region is of high quality (Jones et al. 2003; Pavlenko 
& Jones 2003). In general, flts of our theoretical spectra 
computed with DUSTY model atmospheres to the observed 
spectra look better. In Fig. [5] we show values of min S com- 
puted for a grid of DUSTY and COND model atmospheres. 
It is worth noting: 

• The values of S are smaller for the DUSTY models 
across a wide range of effective temperature and gravity, 
particularly below 2300 K. 

• We get min S at Teii = 2000-2200 K for the COND 
models. This result derived through spectral analysis is con- 
sistent with our estimation of the surface temperature for 
LP 944—20 obtained in Section 3 using radii predicted by 
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Figure 5. Values of min (5) computed for a grid of DUSTY and 
COND model atmospiieres in the "standard" framework. 



substellar evolutionary models, the object's astrometry and 
broad-band photometry. 

• Temperatures in the outermost layers of the DUSTY 
models are higher in comparison with the COND models. 
As s result for the DUSTY models, S appears to decrease 
even for Tcff < 2000K. 

We remark that any spectral analysis is affected by un- 
certainties associated with the adopted structure of model 
atmospheres, molecular opacities, and non-classical effects 
like the presence of a chromosphere, veiling, etc. Addition- 
ally, the fits shown in Fig. 2] do not provide a good match 
to the observed SED of LP 944—20. Therefore, these results 
should be understood as a qualitative confirmation of the 
Teff that we have assumed for LP 944—20. From Fig. g] it 
becomes apparent that additional parameters describing the 
atmosphere of LP 944—20 are needed to obtain a better re- 
production of the observed data. We will further discuss on 
this in Section 5. Next, we will show the results for the high 
resolution "standard" analysis. 
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Figure 6. Fits to the K I and Rb I resonance lines observed 
in the spectra of LP 944— 20 following the "standard" approach. 
The bottom panel shows an enlargement around the Rb I line. 
Theoretical spectrum has solar metallicity. 



4.2 Resonance lines of K I and Rb I 

Resonance doublets of K I and Rb I are recorded simultane- 
ously. This provides the opportunity to study these lines at 
a single epoch. The comparison of CASPEC, SPIRAL and 
UVES spectra shows all these data agree well in spectral 
slope, widths and intensities of lines. The fits to K I and Rb 
I resonance line are shown in Fig. [G] The Rb I line is more 
suitable for analysis, because it is less saturated than the KI 
doublet. We can fit our synthetic spectra to the observed Rb 
I profile taking into account the rotation velocity v sini = 30 
km/s and GFHWM = 0.17 A (see the bottom panel of Fig. 
[B]). Note that other features mainly due to molecular absorp- 
tions of TiO and VO are not so well reproduced in intensity 
by the models. Nevertheless, the Tcff = 2000 K and log g 
= 4.5 model nicely fits the K I and Rb I resonance profiles. 
We note that higher gravities do not provide a better match 
to the observed profiles of the akalis (they predict stronger 
lines), as we will show next. 

Rubidium in LP 944— 20 was also studied by Basri et 
al. (2000). These authors derived Tcff = 2600 K from the 



spectral fitting of the Rb I resonance lines. However, Basri 
et al. (2000) point out that their spectral analysis of the Rb 
I lines indicates higher temperatures than their analysis of 
caesium lines (2400 K). Our computations for 2400 K do 
not provide fits to the SED of LP 944—20 as acceptable as 
those computed for 2000-2200 K (see Fig. [5}. LP 944-20 is 
classified as an M9V field brown dwarf in the literature, thus 
it shares the same spectral classification as LHS 2924, which 
was also included in the sample of Basri et al. (2000) . These 
authors found that LP 944-20 is cooler than LHS 2924 in 
agreement with the more recent results of Dahn et al. (2002) , 
who determine Tcff = 2367 and 2138 K for LHS 2924 and 
LP 944—20, respectively. This value is in much agreement 
with our adopted Toff for the young brown dwarf. We note 
that our synthetic spectra provide reasonable fits to the 
observed profiles of all alkali lines in LP 944—20 using the 
same atmospheric parameters, i.e., without the need of dif- 
ferent values of temperature and gravity for the various lines 
present in the spectra. 
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Figure 7. Fits of DUSTY syntlietic spectra (solar metallicity, 
"standard" approacli) to the observed infrared K I subordinate 
doublet lines. 

4.3 Near-infrared doublet of K I 

We have also fit the NIRSPEC spectrum of LP 944-20, 
which contains the subordinate lines of K I at 1.2432 and 
1.2522 fim. These lines turn out to be gravity-sensitive fea- 
tures as observationally shown by McGovern et al. (2004). 
Their excitation potentials are rather low (~ 1.610 eV), thus 
they can be used for the analysis of stellar spectra. Our fits 
to the observed profiles of K I subordinate doublet using 
the standard modelling procedure are shown in Fig. [T] The 
molecular background at these wavelengths is dominated by 
absorption of H2O ; CrH and FeH contribute to the opacity 
as well. Unfortunately, good quality line lists for CrH and 
FeH do not exist. 

Nonetheless, we get a good fit to the observed pro- 
files of these hues (Fig O using Tcs = 2000 K and log 
g = 4.5. Note that higher gravity models deviate signifi- 
cantly from the observations. From the spectral analysis of 
young M-type objects carried out by Mohanty et al. (2004), 
these authors pointed out that their derived gravities for the 
coolest M dwarfs deviate significantly up to 0.75 dex from 
the isochrone predictions in the sense that measurements are 
smaller. Our gravity derivation for LP 944— 20 also appears 
to be lower than the prediction of log g ~ 5.0 based on the 
evolutionary models by Baraffe et al. (1998) and the likely 
age of the brown dwarf. 

While the overall profiles of the near-infrared K I lines 
are reasonably reproduced by our computations, the cores 
of the lines in the theoretical spectra are deeper than those 
in the observed spectra. These differences in the intensity 
of the line cores of the saturated lines may be explained 
by non-local thermodynamical equilibrium (NLTE) effects 
or chromospheric effects. However, other possible physical 
explanations can be found, especially those related to the 
depletion of refractory elements into condensates in cool at- 
mospheres (see Section 5). 

4.4 Lithium lines and abundances 

The relative strength of atomic lines around the Li resonance 
doublet and TiO bands are shown in Fig. [S] These computa- 



fc. 0.4 



1 1 1 




1 1 1 




V — ''^\ 


A 




TiO Off 






J 1 


TiO on 


1 1 1 -T 1 1 -1 1 1 



0.6690.6695 0.67 0.67050.6710.67150.672 0.67250.6730.67350.674 
Wavelength (micron) 



Figure 8. Theoretical profiles of atomic lines in the 670.8 nm 
region computed both with and without TiO absorption in the 
2000/4.5/0.0 model atmosphere. The strongest line is due to Li 
I. The vacuum scale of wavelengths is used. Synthetic spectra are 
convolved with GFWHM = 1 A and v sini = 30 km/s. Solid and 
dashed lines show spectra computed for log N(Li) = 2.0 and 3.2, 
respectively. 

tions were carried out for two cases to show the dependence 
of the strength of the Li I line on the presence of TiO: (i) in 
the first case we only take into account the absorption for 
atomic lines, and (ii) in the second and more realistic case we 
also accounted for molecular absorption. While in the former 
case, analysis of atomic lines can be carried by equivalent 
width studies, only pseudoequivalent widths (pEWs), i.e., 
equivalent widths (EWs) measured in respect to the local 
psudocontinuum formed by molecular and atomic line back- 
ground, can be obtained from the spectra computed in case 
(ii) (see also Pavlenko 1997). 

Equivalent and pseudoequivalent widths of the Li I res- 
onance doublet differ significantly (see Zapatero Osorio et al. 
2000; Pavlenko 2005). The computed EW of the Li I feature 
at 670.8 nm for the DUSTY model atmosphere 2000/4.5/0 
and log N(Li) = 0.0 is 0.54 and 0.58 A for turbulent ve- 
locities of 2 and 3 km s~^, respectively. For log N(Li) = 
3.25 (i.e., cosmic abundance) we obtain much larger EW of 
13.5 A (this theoretical measurement accounts for the broad 
wings of the doublet ± 40 A from its core). Interestingly, the 
formally measured pEW of the lithium line obtained for log 
N(Li) = 3.25 and the synthetic spectra of case (ii) is consid- 
erably smaller, much in agreement with the observations of 
LP 944—20. Even for case (i) when we only consider atomic 
absorption lines in our computation of synthetic spectra, we 
cannot see the extended but weak wings of the lines due to 
the blending effects by other features. 

From the observations of LP 944-20, the pEW of the 
Li resonance doublet in the UVES spectrum shown in Fig. [5] 
is 0.65 ± 0.05 A. Our pEW measurement is consistent with 
that (pEW = 0.53 ± 0.05 A) of Tinney (1998), suggesting 
little variability in the overall strength of the line. We note 
that these pEWs are very similar to the Li I measurements 
obtained for a Orionis members (3 Myr) of related spec- 
tral types by Zapatero Osorio et al. (2002). Using curves of 
growth measured by these authors for the young M-dwarfs 
in a Orionis (see Fig. 16 in Zapatero Osorio et al. 2002), 
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Figure 9. Fits of the lithium 670.8 nm resonance doublet com- 
puted for DUSTY model atmosphere 2000/4.5/0.0 and different 
lithium abundances in the "standard" framework to observed 
spectra of LP 944— 20. The bottom panel shows an enlargement 
around the lithium atomic line. 



we obtain a lithium abundance of log N(Li) = 3.2 ± 0.3 for 
LP 944—20. Interestingly, the pEW of the Li I hne computed 
for log N(Li) = 3.25 and derived from the synthetic spec- 
trum shown in Fig. [9] (which accounts for TiO absorption) is 
0.65 A, in agreement with the measurement of LP 944—20. 

Nevertheless, synthetic spectra are the most appropri- 
ate tool for a reliable lithium line analysis. Fits of our theo- 
retical spectra computed for the DUSTY model atmosphere 
2000/4.5/0.0 are shown in Fig. (9] the bottom panel of the 
Figure shows the lithium line in larger scale. Theoretical 
spectra were computed for several lithium abundances, from 
depletion by a factor of 10 up to a rich content of log N(Li) 
= 3.5. The best fit is provided by log N(Li)=3.25 ± 0.25 (see 
details in the bottom panel of Fig. [9}. LP 944— 20 seems to 
have preserved its original lithium content. 

It is worth noting that (i) after appropriate broadening 
(see Section 13. 3p the intensity and profiles of the Li I res- 
onance doublet can be well fit, and (ii) the shape of TiO 
bands around 670.8 nm cannot be reproduced with stan- 
dard models. Theoretical spectra predict stronger TiO bands 
than those observed in LP 944—20 (we note that a change 
of 0.5 dex in gravity does not have a major impact on the 



TiO bands in this wavelength range) . In the next section we 
attempt to improve the the poor fit to the TiO bands by 
adjusting our model atmosphere in a manner which seems 
physically justified. 

Lithium subordinate lines at 812.6 and 610.3 nm can 
also be used for Li abundance determination in atmospheres 
of some lithium-rich late-type stars. Moreover, any abun- 
dance derivation becomes more solid if measurements are 
obtained from as many lithium lines and lithium-bearing 
species as possible. We find the Li line at 610.3 nm "sinks" 
under the TiO bands in the LP 944—20 spectrum. Only the 
812.6 nm line can be used for the analysis if reasonably 
good quality high-resolution spectra were available at these 
wavelengths. The CASPEC spectrum of LP 944-20 does not 
comply this requirement. Therefore, in this paper our abun- 
dance analysis will be based on the fits to the observed Li I 
670.8 nm line profile. 



5 ANALYSIS WITH A SEMI-EMPIRICAL 
MODEL ATMOSPHERE 

In the framework of the standard model we have determined 
the "cosmic" lithium abundance [log N(Li)=3.25 ± 0.25] 
in the atmosphere of LP 944— 20. The same result was ob- 
tained from the spectral analysis of the lithium resonance 
doublet profile and the study of the line pEW. However, 
there remain a number of problems. Our fit to the observed 
spectral energy distribution of LP 944—20 is not adequate 
in both the optical and infrared regimes. A comparison with 
high resolution UVES spectra shows significant differences 
in the shapes of computed and observed spectra. Namely, 
the heads of TiO bands look "smoothed" in the observed 
spectra. We find our fits obtained in the framework of the 
standard model cannot be improved by changes of T^s , log 
g , metallicity or microturbulent velocity. 

Generally speaking, these differences reduce the relia- 
bility of our analysis, but on the other hand we know dusty 
effects can affect the temperature structure and spectra of 
ultracool dwarfs (Tsuji et al. 1996; Jones and Tsuji (1997); 
AUard et al. 2000; Burrows et al. 2002b, and references 
therein). Therefore we construct the new semi-empirical 
(SE) model atmosphere. With the new SE model atmo- 
sphere we obtain better fits to the overall spectral energy dis- 
tribution and molecular features around lines of K, Rb and 
Li resonance doublets. Finally, we re-determine the lithium 
abundance with the new model atmosphere. 



5.1 Semi-empirical model atmosphere 

We modified the DUSTY models as follows: 

• To decrease the intensity of TiO bands in the spectral 
range 600-800 nm, we suggest the complete absence of TiO 
above a certain height in the atmosphere of LP 944—20. This 
can have at least two different physical explanations: (i) the 
presence of a hot chromospheric-like region in the outermost 
layers of the dynamical atmosphere, or (ii) the removal of 
Ti atoms from gaseous species by condensation into dust 
particles in the upper atmosphere (see Pavlenko 1998a; Fer- 
guson et al. 2005). In our SE model we will explore the latter 
case, i.e., a depletion of TiO above a particular level in the 
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Figure 10. Structures of two semi-empirical model atmospheres 
from our grid. 



atmosphere of LP 944—20. Our lower boundary of the de- 
pleted TiO is located at the level Tross = 10~*. The choice 
of this point has some influence on the absolute values of 
min S, however, it is not crucial for our results of lithium 
abundance determination. 

• To improve the flts over the wavelength range 0.8-0.9 
and 1.6-2.5 fim , we implement additional quasicontinuum 
opacity (hereafter AqO) provided by dust particles. It is 
worth noting that dust particles must be present in atmo- 
sphere of LP 944—20 due to the low effective temperature 
of the dwarf. Our knowledge about dust opacity in atmo- 
spheres of ultracool dwarfs is rather poor (see, however, 
review by Ferguson et al. 2005). Here we treat AqO as a 
scattering with power-law dependence on wavelength in the 
form ax = ao X (Ao/A)^. We used Ao = 0.77 fim as our ref- 
erence wavelength. Parameter N varies in the range 0-4. A 
value of = corresponds to the hght scattering by sub- 
micron particles ("white scattering"), and A'' = 4 describes 
Rayleigh scattering. 

• In our case, the "dust opacity" originates in the shell- 
like structures laying above the photosphere (clouds). It is 
worth noting that we tried different heights for the location 
of the "dusty shell" . The best results are obtained for levels 
corresponding to Ta ~ 2000 K, coinciding with the models 
of Tsuji (2002) and AUard et al. (2001). These authors re- 
fer to Tcr as a temperature related to the "gas-dust" phase 
transition. 

Our scattering dusty cloud structure was modelled by 
two parameters: thickness of cloud, TAqO , and the location 
of the maximum of opacity, Td . The dusty opacity decreases 
inward and outwards in the atmosphere with respect to the 
Td point (see Fie:.llO[) for physical reasons: pressure and den- 
sity are reduced outwards while temperature increases in- 
wards. These factors naturally decrease the density of the 
dusty particles at lower and larger optical depths. We vary 
TAqO and Td across the ranges 0-1 and 0.001-0.0001, respec- 
tively. Note that different locations of the scattering layer 
correspond in our case to different depths in the atmosphere 
in which T — Tcr- Burrows et al. (2006) have recently pro- 
posed a similar model of the dusty cloud with a flat part of 
the dust distribution located at T ~ 2300 K. 




Wavelength (micron) 

Figure 11. Comparison of fluxes computed with the DUSTY 
2000/4.5/0.0 model atmosphere for two different cases: the dust- 
free atmosphere (no AqO) and the SE model (rAqO = 0.9, Td = 
0.001, k = 25). See also Fig. 2 of Burrows et al. (2006). 

5.2 Fits to the observed spectral energy 
distribution 

We have computed a grid of theoretical spectra adopting 
our SE model atmospheres built with different parameters 
of TAqO , Td , and k (k stands for different layers within 
the atmosphere). A comparison of computed optical and 
near-infrared SEDs for the SE model atmospheres with AqO 
(TAqO =0.9, Td =0.001) and without AqO is shown in Fig. 1111 
As we see from the comparison, our AqO affects the spectral 
regions where fluxes are formed in the layers below the level 
T < TAqO . Strong absorption features formed in the optical 
spectrum are affected rather marginally. In general, the AqO 
makes the computed spectra appear "shallower" than mod- 
els without AqO, i.e., in agreement with the observations of 
LP 944-20. 

Theoretical SEDs computed for SE model atmospheres 
were fit to the observed fluxes of LP 944—20 following the 
procedure described in section 13.51 It is worth noting the 
following: (i) the minimisation procedure of fits to observed 
spectra provides lower S for theoretical spectra computed 
with A'^ = (this suggests that the dusty particles present 
in the modelled dusty cloud structure of LP 944—20 are of 
submicron size), and (ii) visual comparison of Figs. |4] and 
1121 shows that we get a better fit to the observed SED with 
theoretical spectra computed with AqO. Analysis of the re- 
sults of the numerical procedure of minimisation yields S = 
1.4x10"*^ and 3.9x10"'' for the fits shown in Figs. |l(stan- 
dard approach) and [T5] (SE model) respectively. We note 
that while our SE computations decrease the absolute val- 
ues of min S, the general shape of the min S dependence 
on TcfT (like the one shown in Fig. [S| is not significantly 
changed. SE models of Toh = 2000 K provide reasonable 
fits to the observed data of LP 944—20, indicating that our 
results are stable. 

5.3 K I and Rb I resonance lines 

Fits of theoretical spectra computed for SE model atmo- 
sphere to the spectral region of K I and Rb I line are shown 
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Figure 12. Fit of tlie tlieoretical fiuxes computed witii a SE 
model atmospliere (2000/4.0/0.0, TAqO = 0.9, = 0.001, k = 
25) to tlie observed spectral energy distribution of LP 944—20. 
Note the better fit as compared to Fig. |4] 

in Fig. [I3l We get even better fits to K I and Rb I resonance 
lines and to TiO molecular features in terms of intensity in 
comparison to the dusty free case (see Fig. [6]). 



5.4 Near-infrared doublet of K I 

The NIRSPEC data of LP 944-20 is compared to our syn- 
thetic spectrum computed with the SE model atmosphere 
in Fig. 1141 In general, our SE model atmosphere does not 
have a major impact on the atomic and molecular details of 
these wavelengths: the fit of the molecular features around 
K I lines are of the same quality and the profiles of the K 
I lines are reproduced as in the AqO-free case (see Fig. [7]). 
However, we note that the fits to the observed cores of K 
I are improved in the AqO model. And again, we see that 
models of log g —4.5 provide a much better fit to the atomic 
lines than higher gravity atmospheres. 



5.5 Li I resonance lines 

The SPIRAL and UVES spectra of LP 944-20 around Li I 
670.8 nm are compared to SE models in Fig. 1151 These mod- 
els have been computed for different lithium abundances, 
from log N(Li) = 2.0 up to 3.5. The fits to the resonance 
line profile of Li I are of similar quality to those depicted in 
Fig. "standard models"). However, the TiO absorption fit 
and the overall fit of the spectral region across the Li reso- 
nance doublet is significantly improved in the SE models. 

We carry out our lithium abundance determination for 
the SE model atmosphere. As seen from the bottom panel of 
Fig. 1151 the model computed for log N(Li) = 3.25 provides 
the best match to the observations. Interestingly, our SE 
model does not significantly affect the lithium abundance de- 
termination for LP 944—20. Observed and theoretical spec- 
tra have been normalised over a broad range of wavelengths 
across the lithium region because the synthetic data provide 
a reasonable reproduction of the TiO absorption. We have 
adopted an error bar of ±0.25 dex in our lithium abundance 
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Figure 13. Fits of the synthetic fluxes computed with a SE model 
atmosphere (2000/4.5/0.0, TAqO = 0.9, = 0.001, k =25) to the 
LP 944—20 observed spectrum around K I and Rb I lines. The 
bottom panel shown an enlargement around the Rb I line. 
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Figure 14. Fit of computed SE models (2000/4.5/0.0 and 
2000/5.0/0.0, TAqO = 0.9, Td = 0.001, k =25) to the LP 944-20 
observed infrared K I subordinate doublet region. The core of the 
K I lines are better reproduced than in Fig. [7] Note that log g 
= 4.5 provides a better match to the profiles of these "gravity- 
sensitive" atomic lines. 
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Figure 15. Fits of the theoretical fluxes around Li 670.8 
nm resonance doublet computed with SE model atmosphere 
(2000/4.5/0.0, TAqO = 0.9, = 0.001, k = 25) and different 
lithium abundances to the observed spectrum of LP 944—20. The 
bottom panel shows an enlargement around the resonance line. 
Note the better reproduction of the shape and intensity of the 
TiO absorption as compared to Fig. |9] 



determination; this takes into account the quality of the ob- 
served spectra, Tefj changes of ±200 K, and slight modifi- 
cations in the normalisation of the spectra. From our com- 
putations and at the 4 a confidence level we can discard 
lithium depletion factors larger than 10 in the atmosphere 
of LP 944-20. 



6 DISCUSSION 

Our main result has been obtained from the spectral anal- 
ysis of the high-resolution profile of the Li I doublet reso- 
nance fine at 670.8 nm using T^ff = 2000 K {T^s = 2040 K 
is derived for LP 944— 20 from photometric and astrometric 
considerations), log g = 4.5, solar metallicity, and differ- 
ent atmospheric lithium contents from depletion of about a 
factor of 10 up to super-cosmic abundance. We also model 
other alkali lines of K I and Rb I located in other optical and 
near-infrared spectral regions as a test of consistency of our 
results. We note that profiles of the alkali lines obtained on 
different occasions agree quite well in intensities and widths. 
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Figure 16. Location of LP 944— 20 on the lithium depletion 
curves by Baraffe et al. (1998). We have adopted the age de- 
termination of Ribas (2003). 



Some discrepancies may be explained by the signal-to-noise 
differences or variability of dust shell within the atmosphere 
of LP 944-20. 

Firstly, we reproduce the overall spectral distribution 
and profiles of spectral lines of strong absorption lines of 
rubidium and potassium. Second, we find the lithium abun- 
dance in the atmosphere of LP 944-20 by fitting a DUSTY 
model atmosphere to the Li 670.8 nm line profile in the ob- 
served spectrum. We find that the fit to the observed fluxes 
provides a cosmic value of lithium abundance and that anal- 
ysis of the the observed pEW also yields the same log N(Li) 
= 3.25 ± 0.25 dex. 

Tinney (1998) determined a lithium abundance of log 
N(Li)=0.0 ± 0.5 from analysis of the EW of the lithium 
resonance doublet. We can reproduce Tinney's (1998) re- 
sult if we use pEW as EW (see Section [4.41) . as Tinney was 
forced to do by the then available atmospheric models. In- 
deed, only pEW of atomic lines can be measured in respect 
to the background formed by the haze of TiO lines (Pavlenko 
1997, 1997a; Zapatero Osorio et al. 2002). The pEW values 
depend on the strength of the background molecular bands 
of TiO, on the spectral resolution of the data, and on other 
broadening parameters (see Martm et al. (2005) for an ex- 
ample of how rotation affects the equivalent widths of K I 
resonance lines) . Reducing the TiO absorption in the region 
across the 670.8 nm line increases the relative strength of 
the computed Li lines and the formally determined lithium 
abundance decreases. 

It is worth noting that we get some problems with fits 
of theoretical spectra computed for the conventional, i.e., 
self-consistent, model atmosphere to the observed spectral 
energy distribution of LP 944—20. The shape and intensity 
of the TiO bands around the resonance line of lithium are 
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Figure 17. Location of LP 944—20 in tlie H-R diagram and com- 
parison with the evolutionary models of Baraffe et al. (1998). 
Isochrones (solid lines) correspond to the following ages: from 
top to bottom, 100 Myr, 300 Myr, 800 Myr, and 5 Gyr. Evolu- 
tionary tracks (dotted lines) have masses labelled in Jovian units 
(1 Mq ~ 1000 Mj,p). 



not well matched by "standard" models. Generally speaking, 
because of the low surface temperature and activity prop- 
erties of LP 944—20, it might be surprising to obtain good 
fits for conventional models, given that LP 944—20 provides 
substantial evidence of stellar activity phenomena in its at- 
mosphere (Ha and radio emission, strong X-ray flares) . Our 
knowledge of these processes is far from complete even in the 
case of the Sun. Indeed, the presence of these phenomena is 
clear evidence of the vertical stratification of the LP 944—20 
atmosphere. On the other hand, it is plausible that TiO is 
absent in the upper atmosphere of LP 944—20 due to deple- 
tion of Ti atoms into dust particles and/or complete disso- 
ciation of TiO due to higher temperature there. 

To investigate the impact of different factors on our 
main result and in order to provide a better fit to the 
observed spectra of LP 944—20, we adopt a more sophis- 
ticated model atmosphere. Namely, we suggest the pres- 
ence of clouds at the level above the photosphere and com- 
plete depletion of TiO in the outermost layers of the at- 
mosphere of LP 944— 20. We treat the opacity provided by 
the dusty clouds as scattering with power dependence on 
wavelength. From the comparison of observed and computed 
spectra we determine parameters of scattering clouds. The 
temperature-pressure dependence in our model atmosphere 
was taken from the DUSTY model atmosphere of the Lyon 
group, i.e., our SE model is not self-consistent. Even in the 
Sun (e.g., HOLMU, VLC, MAKKLE), the use of SE model 
atmospheres is well acknowledged (see references in Gib- 
son 1970). Semi-empirical model atmospheres are used to fit 
some observed features using theoretical spectra computed 
with modified model atmospheres. 



In our case the appearance of a dusty shell in the SE 
model atmosphere and/or depletion of TiO should change 
the temperature structure of the outermost atmospheric lay- 
ers. With our SE models we obtain considerably better fits 
between observational and synthetic spectra in terms of the 
TiO bands, alkali lines and SED. From fits to the Li reso- 
nance doublet we conclude that LP 944—20 has likely pre- 
served all its originial lithium content. 

We also conclude that our lithium abundance shows 
a rather weak dependence on the input parameters of the 
models, despite the fact that our log N(Li) determination 
is carried out in the frame of local thermodynamical equi- 
librium (LTE). Moreover, all the alkali lines fitted in our 
paper were computed in the frame of LTE. Fortunately, the 
alkali lines considered in this work are rather insensitive to 
chromospheric-like features (GLF) because processes of their 
formation are mainly controlled by photoionisation (Thomas 
1959). By definition, these lines should show a rather weak 
dependence on temperature structure of the outermost lay- 
ers of atmosphere. Indeed, direct NLTE modelling of forma- 
tion of lithium lines in the atmospheres of ultracool dwarfs 
with GLF shows their rather weak response to temperature 
inversions, in contrast to TiO lines (Pavlenko 1998). Lithium 
lines can be affected only in the case of the strongest GLF 
producing additional flux in the blue part of the spectra, 
but that is unlikely in LP 944—20 given its age and level of 
Hq and X-ray activity. 

Finally, comparison of our revised lithium abundance 
for LP 944—20 with the calculations of lithium depletion by 
Baraffe et al. (1998) shows a good agreement for a mass 
below 0.057 Mq and an age of 320±80 Myr (see Fig. [Tg}. 
This result is marginally consistent with the earlier estimate 
of Tinney (1998), who derived a mass in the range 0.056 - 
0.064 Mq. Using the same evolutionary models, we show 
that this age and a mass ^0.057 Mq are also consistent 
with the location of this brown dwarf in the H-R diagram 
(Fie:. ll7p . Thus, our analysis supports the young age of a few 
hundred million years of LP 944—20 as suggested by Ribas 
(2003), and indicates a lower mass than previously reported 
in the literature. 
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